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ABSTRACT 

We use Keck laser guide star adaptive optics imaging and exploit the magnifying 
effects of strong gravitational lensing (the effective resolution is FWHM « 200 pc) to 
investigate the sub-kpc scale of an intermediate-redshift (z = 0.63) massive early-type 
galaxy being lensed by a foreground early-type galaxy; we dub this class of strong 
gravitational lens systems EELs, e.g., early- type/early- type lenses. We find that the 
background source is massive (M* = 10 10 9 M Q ) and compact (r e = 1.1 kpc), and a 
two-component fit is required to model accurately the surface brightness distribution, 
including an extended low-surface-brightness component. This extended component 
may arise from the evolution of higher-redshift 'red nuggets' or may already be in 
place at z ~ 2 but is unobservable due to cosmological surface brightness dimming. 
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1 INTRODUCTION 

The hierarchical formation scenario of galaxy evolution sug- 
gests that galaxies grow via mergers with other galaxies. The 
early-type galaxies (ETGs) that are the products of these 
mergers follow a set of tight scaling relations, including a 
correlation between th e sizes and stella r masses of ETGs 
at lo w redshifts (e.g., IShen et all 120031 ; iHvde fc Bernard] 
2009). However, it has recently been found that ETGs at 
high redshifts (z > 1.5) show substa n tial deviations from 
this relation (e.g.. iDaddi et alj [20051; ITruiillo et al.l l200rj ; 

van Dokkum et alj|2008l ); these 'red 
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nuggets' have effective radii between 2 and 5 times smaller 
for a given stellar mass than is expected from the local r e - 
M» relation. Recent deep spectroscopic observations have 
found that the stellar velocity dispersions of these objects 
are c onsistent with the inf e rred sizes and stel l ar m asses 
(e.g., ICappellari et all 120091 ; Ivan Dokkum et al l [20091) . al- 
though there is some evidence that many galaxies have had 
their sizes systematic ally under-estimated ICappellari et al.1 
l2009r i. Furthermore, iNewman et al.1 (|2010l ) find that the 
sizes of ETGs at slightly lower redshifts than the red nuggets 
(1 < z < 1.6) are at most a factor of two larger than ETGs 
in the local universe, which suggests that most of the size 
evolution of red nuggets - if required - occurs early in the 
lives of these galaxies. 

Hierarchical growth through 'dry' mergers should lead 
to most red nuggets growing substantially in size but not 
as quickly in mass, therefore moving these galaxies towards 
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the local r e -M* relation, although the tight scatter ob- 
served in the r e -M* relation indicates that at most only half 
of th e growth can be due to major mergers (|Nipoti et al.l 
2009). Additionally, up to 10% of the red nuggets will not 
have exp erienced any major m ergers to aid in this 'puff- 
ing up' (|Hopkins et al.l l2009al ) and should thus exist at 
low redshifts. Even if some red nuggets exist locally, most 
will have disappeared, presumably e volving into the cores 
of massive early typ e galaxies fe.g.. IHopkins et al.l [2009b; 
iBezanson et al.l2009l ). Several groups have attempted to find 
the surviving counterparts to these high-z red nuggets at 
low redshifts with v aryin g and sometimes discrepant results. 
IValentinuzzi et all l|2010f ) examined the ETG population of 
the WINGS clusters and found that ~ 20% of these galax- 
ies have properties similar to red nuggets and suggested 
that much of the observed evolu tion may be due to sam- 
ple selection effect s. Conversely, ITruiillo et ail (I2009D and 
iTavlor et ail l|2010l ) searched the SDSS and found that the 
number density of red nuggets must decrease by a factor of 
~ 5000 between z ~ 2 and z « 0.1. This discrepancy makes 
it unclear how many red nugget analogs exist in the local 
universe. 

Most studies have investigated massive red nuggets in 
the distant (z ■> 1.5 ) or local (z « 0.1) universe (but see 
IStockton et al.|[2010l ). This is not due to a lack of interest in 
intermediate re dshifts (0.1 < z < 1); indeed, the rate of evo- 
lution found by iBernardil l|2009l ) suggests that the progeni- 
tors of z = 0.1 BCGs might still be red nuggets at z < 1, and 
IHopkins et al.l (2010) suggest that tracing the evolution of 
the faint 'outer envelope' of ETGs is a strong model discrim- 
inator. The problem with moving to intermediate redshifts 
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has been in identifying compact red galaxies at these red- 
shifts. For example, 0.5 kpc at z = 0.5 corresponds to 0.08"; 
such galaxies may be unresolved even in HST surveys. 

We have implemented a novel approach to find and 
study intermediate-redshift analogs to the high-z red 
nuggets: we use the SDSS spectroscopic database to find 
strong gravitational lens candidates that consist of an 
intermediate-redshift early-type galaxy being lensed by 
a lower-redshift early-type galaxy (we call these early- 
type/early-type lenses, or EELs). These natural telescopes 
provide a unique opportunity to identify and study red 
nuggets at z ~ 0.6, as high mag nification is easier to 
achieve for compa ct sources (e.g., iMarshall et"afl 120071 ; 
iNewton et "ail |201Ch and ETGs are massive. Furthermore, 
the lensing magnification significantly increas es the effec - 
tive depth and resolution of observations (e.g.. iTreul boiO). 
opening the possibility to explore galaxies of unprecedented 
compactness. 

This letter presents the discovery and first analysis of an 
EEL and demonstrates the unique utility of our approach. 
We present the data in Section 2, followed by our lensing and 
surface brightness fitting in Section 3, and conclude with a 
discussion in Section 4. All magnitudes are in the AB system, 
stellar masses are calculated assum ing a Salpeter initial mass 
function fe.g-. lAuger et alj|2010al ). and we use a flat ACDM 
cosmology with fi m = 0.3 and h — 0.7. 



2 OBSERVATIONS AND DATA REDUCTION 

We have used the SDSS spectroscopic database and archive 
to identify spectra that appear to be composed of two early- 
type galaxy spectra at different redshifts and are therefore 
candidate strong gravitational lenses; this selection is sim- 
ilar to that emplo yed by the Sloan Lens ACS (SLACS; 
iBolton et all [2006) survey, although the SLACS survey 
looked for emission lines from the background objects and 
therefore selected late-type sources. The selection method 
essentially guarantees that the background galaxy is intrin- 
sically luminous and therefore massive, as the spectroscopic 
continuum must be clearly visible in order to distinguish 
both the foreground and background redshifts. 

We have observed 7 EEL candidates to date using 
NIRC2 with the laser guide star adaptive optics (LGS- 
AO) system on Keck II and have successfully confirmed 
that all of the observed can didates are strong gravitational 
lenses (| Auger et alj|2010bh . Here we focus on the system 
SDSSJ1347-0101 with lens deflector redshift z = 0.39 and 
lens source redshift z = 0.63, which was observed on 21 
April 2010 UT in clear skies using a tip-tilt star with mag- 
nitude r — 16.8 offset by 56"from the lens system. We used 
the NIRC2 wide camera (with pixel scale / .'04 pixel -1 and 
O'.'ll FWHM resolution during our observation) with the 
Kp filter and obtained 30 exposures of length 40s for a total 
on-source time of 1200s. Each exposure also includes a star 
located 17"from the lens which we use as the PSF star. 

The data were reduced by flatfielding with a sky flat 
created from the median of the science exposures with all 
sources masked, and a median background was then sub- 
tracted from each image. These images were resampled to 
a distortion-corrected frame and the 30 images were then 
registered using cross-correlation and pixel-based techniques 



which find sub-pixel offsets between the images. These off- 
sets and the distortion model were then used to drizzle 
the images onto a common output frame. The registration 
worked very well for the PSF star but we found that the dis- 
tortion model left residual offsets of up to 0.7 pixels at the 
lens location, and we therefore re-registered the images with 
the PSF star masked to obtain an un-smeared final image 
of the lens system (Figure [TJ . 



3 MASS AND SURFACE BRIGHTNESS 
MODELLING 

We model the lens surface m ass distribution as a singular 
isothermal ellipsoid (SIE; e.g.. lKormann et al.lll994[ ) and al- 
low for external shear. The light of the lensing galaxy is 
modelled as an elliptical Sersic profile with Sersic index n 
free to vary between 0.5 and 6, but we do not assume that 
the mass and the light of the lens are aligned. The back- 
ground source was also initially modelled as a Sersic profile 
with Sersic index between 0.5 and 6. We convolve the mod- 
els with an empirical PSF that was observed simultaneously 
with the lensing system and compare the convolved models 
with the data. 

The best-fit model is obtained by using a Levenburg- 
Marquardt optimisation that minimises the x 2 of the differ- 
ence between the data and the model. The data are generally 
well-described by the model, but we find that using a single 
Sersic component for the source results in a ring of resid- 
ual flux left in the image (Figure [TJ ; the amplitude of these 
residuals is approximately 1.5 to 3 times the noise level. We 
also find that the data do not provide precise constraint on n 
for the source (they only constrain n to be large, i.e., n > 3) 
so we fix the source to have an elliptical de Vaucouleurs pro- 
file (n = 4). The residual flux is clearly not associated with 
the foreground deflector, as adding a second Sersic compo- 
nent does not improve the fit. However, adding a second 
component to the source surface brightness profile improves 
the fit dramatically, as shown in Figure [T] and the per-pixel 
significance of the residuals around the ring is below the 
noise. The best-fit parameters for the lensing mass distribu- 
tion and the light profiles for the source for both the one- 
and two-component models are given in Table [1] 

We find that the foreground galaxy is well-described by 
a single Sersic component for the light a SIE for the mass dis- 
tribution. The SIE Einstein radius is CK'42, which implies a 
lens velocity dispersion of 210±15 km s -1 . This is consistent 
with the stellar velocity dispersion of 203±42 km s _1 that 
we ob tain by fitting velo c ity-di spersion-broadened single- 
burst iBruzual fc Charlotl (120031 ) templates to the SDSS 
spectrum; Figure [5] shows the result of this fit, which we 
use to decompose the spectrum into lens and source com- 
pone nts. The observed l e ns Kp magnitude is 17.7 and we 
use a lBruzual fc Cha riot (2003) instantaneous-burst galaxy 
template with formation redshift z — 3 and solar metal- 
licity to determine rest-frame quantities. We find a K- 
band rest-frame magnitude of 18.4 and a luminosity of 
= 10 Lk.q, with a corresponding stellar mass 
M*,i cns = 10 11,6 M . We adopt error estimates of 0.1 dex on 
the luminosity and 0.2 dex on the stellar mass for both the 
lens and source. 

The source has a total Kp magnitude of 20.3 (20.6 
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Figure 1. Keck LGS-AO imaging of the EEL SDSSJ1347-0101 (left). The image is well modelled using a singular isothermal ellipsoid 
mass distribution for the lensing mass, a single Sersic component for the lensing galaxy, and a two-component de Vaucouleurs and Sersic 
profile for the background source (right), while a single-component fit to the source leaves a ring of residual flux (centre); the scaling of 
the two residual plots is set to saturate at 3 times the background noise level. The Kp-band imaging reveals an observed source magnitude 
of 17.6 but this is corrected to 20.3 when the magnification from lensing is taken into account; the inferred stellar mass of the source is 
10 10,9 Mq and the half-light radius of the two-component surface brightness model is 0'.'16 or 1.1 kpc at the redshift of the source. 



Table 1. Best-fit parameters for the lens mass model and source light distribution of SDSSJ1347-0101. Columns: (1) Description of the 
source surface-brightness model; (2) lensing Einstein radius in arcsec; (3) lens model velocity dispersion in km s _1 ; (4) lens mass axis 
ratio b/a; (5) lens mass position angle in degrees; (6) external shear; (7) external shear position angle in degrees; (8) source magnitude; 
(9) source effective radius in arcsec; (10) source Sersic index. The errors on mxp are 0.15 mags and the errors on r e are 10%. 
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Figure 2. Smoothed SDSS spectrum of SDSSJ1347-0101 (black) 
and model template spectra for the lens (blue) and source (red). 
Ca K &; H lines from the background source are clearly visi- 
ble in the smoothed spectrum, and we use the lens/source de- 
composition to determine a stellar velocity dispersion cgrjss = 
203 ± 42 km s~ 1 for the lens. Furthermore, we find that the stel- 
lar masses required to normalise the template spectra to match 
the observed spectrum are consistent with the masses derived 
from the Keck AO imaging (see Section 3). 



if a single-component is fit) and, again assuming a for- 
mation redshift z = 3 and solar metallicity, this implies 
that the rest-frame K-band magnitude is 21.0, the lumi- 



nosity is LK,src = 10 Lk,q, and the stellar mass is 
M„, src = 10 10 - 9 M . The lensed (i.e., observed) source mag- 
nitude is 17.6, and we therefore find a factor of 12 magnifi- 
cation due to lensing. The background source size is 0.5 kpc 
if a single de Vaucouleurs profile is assumed or 1.1 kpc if a 
two-component de Vaucouleurs and Sersic profile is used. 



The SDSS spectrum allows us to verify that our stellar 
m asses from the Keck AO-LG S photometry are robust. We 
fit lBruzual~fe Chariot] (|2003l ) templates to the foreground 
and background components of the spectrum to determine 
the appropriate normalisation of the templates and therefore 
the stellar mass of the flux observed by the SDSS spectro- 
scopic fibre. The lens (foreground) stellar mass determined 
from the SDSS spectrum is 10 11 ' 1 Mq which implies a 0.5 
dex loss of flux due to the fibre aperture; this is consis- 
tent with what we find if we convolve our model with the 
SDSS point spread function (which has FWHM = 2" 2 for 
SDSSJ1347-0101) and integrate the flux within the SDSS fi- 
bre aperture. We perform the same procedure for the lensed 
(background) source and find that, when accounting for 
magnification and fibre losses, the source stellar mass in- 
ferred from the Keck LGS-AO photometry and the SDSS 
spectrum agree to within 0.2 dex, consistent with our as- 
sumed errors on the stellar mass. 
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Figure 3. The size -mass relation f or hig h-redshift ETGs; black 
diamonds are from iNewman et al.l ll2010h . grey circ les are from 
ICimatti et alJ |2008f). and light-g r ey squ ares are from lKriek et all 
1120081 ) and Ivan Dokkum et al.1 j2008h . SDSSJ1347-0101 (red 
points; the square is for a single-component de Vaucouleurs fit, 
the hexagon is for the two-component fit, and the triangle is the 
de Vaucouleurs component of the two-component fit) has proper- 
ties consistent with the high-redshift objects but is a significant 
outlier from the local-universe size-m ass relation (black line w ith 
1-cT scatter indicated by dotted lines; iHvde fc Bernardill2009r) . 
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Figure 4. The observed surface density profile of the background 
source of SDSSJ1347-0101 (dotted grey line) and best-fit one- 
component (green) and two-component (purple) models plotted 
from the HWHM of the effective resolution to the surface bright- 
ness limit of our observations. We also plo t the surface br i ghtness 
profiles for massive Virgo ellipticals from [Hopkins et al. (2009b) 
rescaled to a Salpeter I MF (blue lines) and for t he z ~ 2 red 
nuggets (red lines) from Ivan Dokkum et al. | ll2008l l. The single- 
component model is unable to fit the low-surface-density tail and 
ovcrfits the central region to compensate, as is seen in Figure 1. 



4 DISCUSSION AND CONCLUSIONS 

The background source of SDSSJ1347-0101 is a clear out- 
lier from the size-mass relationship of local ETGs and has 
properties that are similar to high-redshift ETGs (Figure [3]). 
However, the inferred size varies by a factor of two depend- 
ing on whether a one- or two-component surface brightness 
model is fit to the galaxy; our imaging data clearly indi- 
cate that the two-component model is required to describe 
the light distribution of SDSSJ1347-0101. Larger samples of 
intermediate-redshift objects are need to make any strong 
claims, but we note that the extended 'wing' of flux from 
the background source of SDSSJ1347-0101 may help alle- 
viate the tension between the very compact high-redshift 
galax ies and the apparen t lack of local-unive rse counterparts 
(e.g.. lTruiillo et al.ll2009l ; iTavlor et alj|20ich in two ways. 

First, single-component fits may be poor models for 
the t rue surface brightne ss distribution of compact ETGs 
(e.g.. lStockton et al.l20l6T ). and observations of high-redshift 
galaxies may lack the sensitivity required to characterise 
low-surface-brightness extended components of red nuggets. 
This is unlikely to fully resolve the issue, however, since a 
factor of two size difference is not enough to move all of the 
high-redshift objects to the local size-mass relation. We also 
note that deep HST imaging of the high-redshift galaxies 
does not indicate the presence of a significant second co m- 
ponent for most red nuggets (e.g.. iDamianov et al.ll2009l ). 

Alternatively, the low-surface-brightness components of 
intermediate-redshift red nuggets mi ght be due to evolution 
of the higher-redshift red nuggets. iHopkins et alJ |2009b) 
show that massive ETGs in the local universe have ex- 
tended low surface brightness components but the central 
surface mass density of these galaxies is similar to the cen- 
tral surface mass density of the red nuggets, indicating that 



the red nuggets may 'grow' extended wings of low-surface- 
brightness material. We find that SDSSJ1347-0101 is consis- 
tent with this scenario and has a peak (PSF-deconvolved) 
surface mass density of lO ll :2 M0kpc _2 , similar to both low- 
and high-redshift massive ETGs (Figure [4}. 

IStockton etTaO l|20ld ) used 9-band photometry from 
SDSS and UKIDSS to identify compact objects that have 
colours consistent with ETGs at redshifts 0.4 < z < 0.6. 
They found 14 compact massive galaxy candidates in an area 
of ~200 square degrees and followed up 11 the them; 7 of 
these 11 were spectroscopically found to be galaxies and two 
were confirmed to be compact (r e ~ 1 kpc) using LGS-AO 
imaging with NIRC2. These intermediate-redshift, massive 
compact galaxies also require low-surface-brightness compo- 
nents to adequately describe their surface brightness distri- 
butions, and this may indicate that lower-redshift compact 
massive galaxies require two-component models in general. 

A larger number of intermediate-redshift red nuggets 
must be investigated to determine quantitatively their den- 
sity, test the degree to which different physical mechanisms 
affect the evolution of the high-redshift red nuggets, and 
characterise the scatter in the si ze-mass relationship of these 
objects (e.g.. iNipoti et al1l2009l ). Depending on the slope of 
the size-luminosity function, lensed sources may be more 
compact than non-lensed counterparts due to the increased 
cross- section for lensing of high-surface-brightness sources 
(e.g.. iNewton et al] |2010). Our Keck adaptive-optics imag- 
ing has been extremely successful at confirming the lensing 
nature of EEL candidates discovered in the SDSS spectro- 
scopic database. However, there are several potential draw- 
backs to this approach, all of which are surmountable. 

The first is the non-trivial selection function which 
makes calculating precise number densities, for example, 
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difficult. iDobler et alj (|2008l ) have characterised the selec- 
tion function for the SLACS survey, and a similar approach 
will yield the appropriate volume in which we search for 
EELs. We can nevertheless obtain an order of magnitude 
estimate of the number density of EEL sources by mak- 
ing some approximations. We find that roughly one out 
of every 1000 SDSS spectra (which have circular apertures 
with radius l'.'5) contain a high-quality EEL candidate. If 
we assume that we are uniformly sensitive to sources be- 
tween 0.3 < z < 0.7 (this is an approximation; we ne- 
glect magnification bias and the actual selection is biased 
towards the higher redshift end of this interval due to the 
lensing cross section), this yields an EEL source density 
of 1.3 x 10 -3 Mpc -3 . This is a factor of a few smaller 
than the density of massive ETGs at these redshifts (e.g., 
iBundv et al]|2006h . indicating that EELs sources are not an 
extreme subset of the overall population For comparison, the 
number density of EELs s ources is slightly hi gher than high- 
redshift red nugg ets (e.g.. lTavlor et alj|20ich and is several 
orders of magnitude larger tha n that of intermediate- redshift 
red nuggets inferred from the IStockton et~aH (|2010r i search 
(they find one object per 100 square degrees at 0.4 < z < 0.6 
while we find one EEL in every 2 square arcminutes at 
0.3 < z < 0.7). This difference is due in part to our higher 
sensitivity (and therefore lower stellar mas s limit), but it 
also in dicates that the selection function of IStockton et all 
(|201Gl ) is much more restrictive and suggests that the EEL 
source population may be more representative of the overall 
population of compact red galaxies at intermediate redshifts. 

Two other drawbacks of our approach are related to the 
use of AO imaging to confirm and analyse the sources. The 
requirement of having a bright tip-tilt star within one ar- 
cminute of the lens substantially limits the ability to follow- 
up targets and therefore decreases the possible sample size. 
Additionally, accurately modelling the sources requires a 
good model for the AO PSF; in principle this can be ob- 
tained by observing stars off-axis from the science observa- 
tion, although we have found that results are significantly 
more robust when the PSF star is observed simultaneously 
with the lens, as in the present case of SDSSJ1347-0101. We 
note, however, that space-based imaging of EEL candidates 
alleviates both of these concerns due to the stable PSF from 
space and the ability to observe any part of the sky. 

Finally, we note that the strong lensing magnification 
increases the observed flux from the lensed source and there- 
fore makes it easier to obtain spectroscopic follow-up of the 
background source (although the presence of the foreground 
lensing galaxy complicates the interpretation of the spectra). 
Our program therefore provides an efficient means of explor- 
ing the size-mass and mass-velocity dispersion relations of 
massive compact galaxies at intermediate redshifts, and we 
may be able to use spectroscopy to explore the stellar pop- 
ulations of the sources in detail to investigate evidence of 
a younger population of stars associated with the extended 
low-surface-brightness profile. 
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